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Gibbons et al. [Phys. Rev. A 70, 062101 (2004)] have recently defined a class ol discrete Wigner 
functions W to represent quantum states in a finite Hilbert space dimension d. I characterize the set 
Cd of states having non-negative W simultaneously in all definitions of W in this class. For d < 5 
I show Cd is the convex hull of stabilizer states. This supports the conjecture that negativity of W 
is necessary for exponential speedup in pure-state quantum computation. 
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I. INTRODUCTION 

Continuous-variable quantum systems can be repre- 
sented in phase space using various quasi-probability dis- 
tributions, notably the Wigner function W{q,p) 0, 
This real-valued function plays some of the roles of the 
classical Liouville density, for example allowing us to 
calculate some system properties through phase-space 
integrals weighted by W{q,p). Despite these similari- 
ties, W{qo,po) cannot be interpreted as the probability 
of simultaneously measuring observables p and q with 
eigenvalues pq and qo: such dispersion-free values for 
non-commuting observables are not allowed in quan- 
tum mechanics. In fact, W{q,p) can even be negative 
in some phase-space regions, something that obviously 
could not happen were W{q,p) a true probability distri- 
bution (hence the term gwasz-probability) . 

The connection between negativity of W{q, p) and non- 
classicality has not been completely fleshed out, partly 
due to different subjective views on what qualifies as 
'non-classical behavior' (see |EB)- Negativity of W^(g,p) 
has been linked to non- locality but the relation is 
not straightforward. For example, the original Einstein- 
Podolsky-Rosen state can display non-locality despite 
having positive W{q,p) in all phase space 

Buot ^\ and Hannay and Berry T^l seem to have been 
among the first to propose analogues of the Wigner func- 
tion for finite-dimensional Hilbert spaces. Their finds 
were rediscovered later by Cohen and Scully 
Feynman J^] , who defined a discrete Wigner function W 
for the case of a single qubit. This work was developed 
by Wootters [l3| and Galetti and de Toledo Piza [l3| . 
who introduced a Wigner function for prime-dimensional 
Hilbert spaces. There followed other definitions valid for 
dimension d which is odd , even |l6j , power-of- prime 
[Ullil, or arbitrary O Hfl HH. These have been re- 
cently used to visualize and get insights on teleportation 
123, , quantum algorithms i^ij , and decoherence 
|2^. A recent review of phase-space methods for finite- 
dimensional systems is given in |25| . 

In this paper I investigate the relation between neg- 
ativity of discrete Wigner functions and quantum com- 
putational speedup. I will focus on the class of Wigner 



functions defined by Wootters 0| and Gibbons et al. 
|l8j | for power-of-prime dimensions. Wigner functions in 
this class are defined by associating lines in a discrete 
phase space to projectors belonging to a fixed set of mu- 
tually unbiased bases. We will see that the set of states 
displaying negativity of W depends on a number of ar- 
bitrary choices required to pick a particular definition 
of W from the class. I eliminate this arbitrariness by 
characterizing the set Cd of states having non-negative 
W simultaneously in all definitions in the class, for d- 
dimensional Hilbert space. 

For a single qubit, I show that the set C2 consists of 
states which fail to provide quantum computational ad- 
vantage in a model recently proposed by Bravyi and Ki- 
taev [23. For dimensions 2 < d < 5, I show that the 
set Cd is the convex hull of a set of stabilizer states, i.e. 
simultaneous eigenstates of generalized Pauli operators 
[27L l23 |. This is interesting, as quantum computation 
which is restricted to stabilizer states and gates from the 
Clifford group can be simulated efficiently on a classical 
computer [23. If the result holds for arbitrary power- 
of-prime dimensions (as I conjecture), then pure states 
in Cd would be 'classical' in the sense of having an effi- 
cient description using the stabilizer formalism. This sug- 
gests that states with negative W may be necessary for 
exponential quantum computational speedup with pure 
states. 



II. DISCRETE WIGNER FUNCTIONS 



In this section we review the discrete Wigner func- 
tions introduced by Wootters and Galetti and de 
Toledo Piza 0| for prime dimensions, and elaborated 
on recently by Wootters fl^ and Gibbons, Hoffman and 
Wootters UJ for prime-power dimensions. We start by 
defining discrete analogues of phase space and its parti- 
tions into parallel 'lines' (i.e., striations). Then we review 
how to define a class of discrete Wigner functions W by 
associating lines with projectors onto basis vectors from 
a set of mutually unbiased bases (MUB's). 
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A. Phase space and striations 

The discrete analogue ol phase space in a d- 
dimensional Hilbert space is a d x d real array. Unlike 
the continuous phase space, in this discrete setting we do 
not have geometrical lines of points. Instead, a 'line' is 
defined as a set of d points in discrete phase space. 

It is clear that our discrete phase space can then be 
partitioned in multiple ways into collections of parallel 
lines (i.e. disjoint sets of d phase-space points). Following 
Wootters ,17.] . we call each such partition a striation. 

In ^3 a procedure for building (d+l) striations of a dx 
d phase-space array was outlined, resulting in striations 
with the following three useful properties: 

i Given any two points, exactly one line contains both 

points; 

ii given a point a and a line A not containing a, there is 

exactly one line parallel to A that contains a; 

iii two non-parallel lines intersect at exactly one point. 

The construction of these (d -I- 1) special striations in- 
volves labeling the discrete phase space with elements of 
finite fields, and defining lines using natural properties 
of the field (for details, see These striations play a 

central role in the definition of the discrete Wigner func- 
tion W, as we will see below. 



B. Mutually unbiased bases 

Wootters' definition of discrete Wigner functions 
makes use of a special set of d-t-1 bases for a d-dimensional 
Hilbert space. Consider two different orthonormal bases 
Bi and B2: 

Bi ^ {\ai.i) , |ai,2) , \{ai^i\ai,j)\^ = Sij, (1) 

B2 = {|a2,i) ' l"2,2) , \a2.d)}, \{a2,i\'^2.j)\^ = (2) 

These two bases Bi and B2 are said to be mutually un- 
biased, or mutually conjugate, if 



\{atj\ak,i)\'^ = i if i ^ fc. 



(3) 



Wootters and Fields showed that one can define (d + 
1) such mutually unbiased bases (MUB's) for power-of- 
prime dimension d [s^- Note that this is exactly the 
number of striations one can find with properties i-iii 
above; we will use this now to define a class of discrete 
Wigner functions W. 



C. Defining a class of discrete Wigner functions 

We now have the ingredients to define a class of discrete 
Wigner functions: a set of (d -|- 1) mutually unbiased 
bases {Bi, B2, Bd+i}; and a set of (d -I- 1) striations 



{Si, S2, Sd+i} of our d X d phase space into d parallel 
lines (i.e. disjoint sets) of d points each. To define a 
discrete Wigner function we need to choose two one-to- 
one maps: 

• each basis set Bi is associated with one striation 
Si] and 

• each basis vector \aij) is associated with a line A^j- 
(the jth line of the ith striation) . 

With these associations, the Wigner function W is 
uniquely defined if we demand that 



Tr{ 



(4) 



i.e. we want the sum of the Wigner function elements 
corresponding to each line to equal the probability of 
projecting onto the basis vector associated with that line. 

Note that there are multiple ways of making these asso- 
ciations. In general, this will lead to different definitions 
of W using the same fixed set of MUB's. The procedure 
outlined above then leads not to a single definition of W, 
but to a class of Wigner functions instead. 



III. WIGNER FUNCTION FOR A QUBIT 

Let us illustrate Wootters' definition of a discrete 
Wigner function with the simplest case, a qubit. The 
discrete Wigner function W is defined on a 2 x 2 array: 







1^2,1 


1^2,2 



(5) 



There are three striations of this phase space with 
properties and iii as required. Below I list them, 
using numbers j = 1, 2 to arbitrarily label the two lines 



Xij in each striation Si. 



Si : 



1 


1 


, 5*2 . 


1 


2 




1 


2 
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1 


2 


2 


1 



(6) 



Now we need to define a set of 3 mutually unbiased 
bases for a qubit. These can be conveniently chosen as 
the eigenstates of the three Pauli operators ax , and 
az- Let us label these basis vectors \aij), where i S 
{1, 2, 3} indexes the MUB {i — 1,2,3 respectively for the 
operators (Jx,(7y,<Jz), and j € {1,2} indexes the basis 
vector in each MUB (j = 1,2 indicates eigenstates with 
eigenvalues respectively equal to -1-1,-1). 

We can now define W by imposing condition Q, i.e. 
we want the sum of Wij in each line A to be the prob- 
ability pij of projecting the state onto the basis vector 
\a,. 



Pi,j = Tr{\aij) (aijl/5) 



E 



(7) 
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Using the three striations we defined, these conditions 
can be exphcitly stated, in terms of the probabihties pi j: 



Wi.2 



Wl^2 
W2^2 
W2S 
W2,2 
W2,2 
W2.I 



Pi, 2, 
P2,l, 
P2,2, 
P3,l, 
P3,2- 



(8) 
(9) 
(10) 

(11) 
(12) 

(13) 



The equations above define uniquely the Wigner func- 



tion W in terms of the probabihties Pi^j 






W^i.i = \ (Pi.i ~ 


1- P2,l H 


^P3,l - 


1), 


(14) 


Wi^2 = ^ (pi,i - 


\- P2,2 H 


hP3,2 - 


1), 


(15) 


W^2,l = ^ (Pl,2 - 


|-P2,1 H 


KP3,2 - 


1), 


(16) 


W^2,2 = ^ (Pl,2 - 


f P2,2 - 


f P3,l - 


1). 


(17) 



Here we should keep in mind that not all probabilities 
Pij are independent, as the sum over any MUB must add 
up to 1 (i.e. yi,J2jP^,j = !)■ 

From the definitions p4l) - (|17|l we can immediately read 
the conditions for non-negativity of the Wigner function: 



Pl.l +P2,1 +P3,1 > 1 

Pl.l +P2,2 +P3,2 > 1 

Pl,2 +P2,1 +P3,2 > 1 

Pi, 2 +P2,2 +P3,1 > 1 



(18) 
(19) 
(20) 
(21) 



Because of the conditions "YlijPij — 1; there are only 
three free variables in the inequalities above. We can 
thus pick one pi j for each MUB i (say, Pi.i); repre- 
sent states by points in this three-dimensional probability 
space. This representation is equivalent to the more fa- 
miliar Bloch sphere, with the advantage of generalizing 
easily to power-of-prime dimension d. In that case (to be 
discussed in sectional, there are a total of {(P — 1) inde- 
pendent probabilities pij describing any mixed quantum 
state. Note that the probabilities pij combine linearly for 
convex combinations (probabilistic mixtures) of quantum 
states, a consequence of the linearity of the trace in eq. 

O- 

The inequalities (|18|I - H21|I for non-negativity of W are 
satisfied by state p= (pi,i,P2,i,P3,i) if and only if p lies 
inside the tetrahedron 

Ti = convex hull of {(0, 0, 1), (0, 1, 0), (1, 0, 0), (1, 1, 1)} 

(22) 

(see figure ^1. 

By contrast, there exist quantum states whose W is 
negative in some phase-space points, i.e. lying outside of 
tetrahedron Ti. This is illustrated in figure ^ where the 
set of points corresponding to quantum states is the ball 
of radius r = 1/2 and with center (i, i, i). 




0.6 
0.4 n 

0.2 Pi J 



FIG. 1: Tetrahedron in p = (pii,P2i,P3i) space representing 
states with non-negative Wigner function W , for the defini- 
tion of W given by eqs. p4^ - H17j . The ball represents one- 
qubit quantum states. 



Let me draw attention to two features of the set of 
states with non-negative W . First, there exist points p £ 
Ti which do not correspond to any one-qubit quantum 
state (e.g. p = (1, 1, 1)). We must also keep in mind that 
there were arbitrary choices involved in the particular 
definition of W that we picked. We have chosen arbitrary 
one-to-one maps between MUB's and striations, and also 
between lines in a striation and basis vectors. In the next 
section we take these two points into account, defining a 
set which consists solely of quantum states, and for which 
the Wigner function is non-negative for all definitions of 
W. 



IV. ONE-QUBIT STATES WITH 
NON-NEGATIVE WIGNER FUNCTIONS 

In this section I use negativity of W to define the set 
Cd of states in d-dimensional Hilbert space having non- 
negative W for all definitions of W based on a fixed set 
of MUB's. I will then argue that single-qubit states in C2 
behave 'classically' in a concrete computational sense. 

Having fixed a set of 3 MUB's for a single qubit, 
Wigner functions W can be defined in a number of ways, 
corresponding to the 3!(2!)'^ — 48 different associations 
between lines and basis vectors, and between striations 
and MUB's. For each of these 48 Wigner function def- 
initions we can do as above, and find the set of points 
in p-space for which W is non-negative. A simple cal- 
culation shows that depending on the definition, this 
set is either the original tetrahedron Ti or the tetrahe- 
dron whose vertices are the other four p-cube vertices: 
T2 = convex hull of {(1, 1, 0), (1, 0, 1), (0, 1, 1), (0, 0, 0)}. 

Now I define the set Cd of states in d-dimensional 
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FIG. 2: Octahedron in p = {pii , P21 , Psi) space representing 
states with non-negative Wigner functions W for all defini- 
tions of W using a fixed set of mutually unbiased bases. 



Hilbert space which I will argue behave 'classically' in 
a computational sense: 

Definition: The set Cd is defined as the states in a d- 
dimensional Hilbert space whose Wigner function 
W is non-negative in all phase space points and for 
all definitions of W using a fixed set of mutually 
unbiased bases. 

fn the remainder of this article I will characterize the set 
Cd for some small dimensions d, and discuss the limita- 
tions of doing quantum computation solely with states in 
Cd. 

For a qubit, the set C2 is given in p-space as the inter- 
section of the two tetrahedra Ti and T2 presented above. 
This corresponds to an octahedron inscribed inside the 
ball of quantum states: 

C2 = convex hull of {(1, i i), (0, ^ i), 1, i), 

(^,0,i),(i i i 0)} (23) 

(see figure 0). We see that set C2 can be characterized 
in a simpler way as the set of states which are convex 
combinations of the six basis vectors of the three chosen 
MUB's: 

p e C2 ^ p ^^qi,j\aij) {a^j] , ^qi,j=l. (24) 

There are at least two motivations for considering the 
set C2, as opposed to states with non- negative in a 
single definition (the sets Ti or T2). The first one is the 
realization that a priori there is no preferred definition of 



W from the full class of definitions, and any concept of 
'classical states' based on non-negativity of W should be 
definition- independent. The second motivation only be- 
came apparent after we calculated C2' unlike the sets Ti 
and T2, all states in C2 are physical states, i.e. obtainable 
from measurements on a single qubit. 

In the next section we will review other properties of 
C2 from a computational perspective. 

A. States in C2 and quantum computing 

In the last section I defined the set C2 of single-qubit 
quantum states which have non-negative W in all defini- 
tions of W using a fixed set of MUB's. Let us now review 
an argument for the 'classicality' of the set C2, from a 
computational point of view. 

Recently Bravyi and Kitaev j2^ proposed the following 
model of computation. Imagine that for some reason 
there are a few quantum computational operations which 
we can perform perfectly - let us denote these by Oideai ■ 
In addition to those, some operations in set O faulty can 
only be performed imperfectly. Now let us consider a 
particular choice for Oideai- 

• Prepare a qubit in state |0) (i.e., an eigenstate of 
the cTz operator) ; 

• Apply unitary operators from the Clifford group 
(such as the Hadamard and CNOT gates); 

• Measure an eigenvalue of a Pauli operator (dx , d'y 
or ffz) on any qubit. 

The Gottesman-Knill theorem states that the oper- 
ations in Oideai above can only create a restricted set 
of states known as stabilizer states, i.e. simultaneous 
eigenstates of the Pauli group of operators. Moreover, 
such operations do not allow for universal quantum com- 
putation, and can be efficiently simulated on a classical 
computer. 

In addition to these perfect operations, Bravyi and Ki- 
taev proposed a set O faulty with a single extra imperfect 
operation: 

• Prepare an auxiliary qubit in a mixed state p. 

In this model of computation, it is easy to see that 
auxiliary qubits in states p G C2 cannot be used to per- 
form universal quantum computation |2||. States in C2 
are linear convex combinations of the six eigenstates of 
single Pauli operators. As such, they can be prepared ef- 
ficiently from operations in Oideai , namely, Clifford gates 
to obtain any Pauli eigenstate deterministically from the 
initial cTz eigenstate, together with classical coin tosses to 
prepare them with the appropriate weights given by the 
desired convex combination. These are operations that 
can be efficiently simulated on a classical probabilistic 
computer. This justifies my 'classicality' claim for states 
in C2. 
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Interestingly, in '2^ it was shown that some 'non- 
classical' states p ^ C2 can be used to attain universal 
quantum computation in this model. The basic idea is 
to use auxiliary pure states outside of C2 to implement 
gates outside of the Clifford group, using operations in 
Oideai only. A single generic non-Clifford gate together 
with the set of Clifford gates allows for universal quan- 
tum computation. This procedure works also for a large 
class of mixed states p ^ C2, through a distillation of 
pure non-stabilizer pure states from p^^ using Clifford 
operations only [26| . 

This enables us to identify non-stabilizer states as a 
resource which can be tapped to implement non-Clifford 
gates and achieve universal quantum computation. This 
idea was first suggested by Shor 'sT] , and has since been 
elaborated on by other authors (3^ ,3^ i35|, |36|] . 



V. HIGHER DIMENSIONS 

We can follow Gibbons et al. and define a discrete 
Wigner function W whenever the Hilbert space dimen- 
sion d is a power of prime. In this section I use the def- 
inition of Cd to characterize this set for states in small 
Hilbert space dimensions d. 

In d dimensions, the probability-space point p describ- 
ing each state has [d^ — I) components. Requiring non- 
negativity of W for all definitions will correspond to a set 
of inequalities for p, each delimiting a half-space where 
W is non-negative at a particular phase-space point, and 
for a particular definition of Wigner function. States in 
the set Cd are those which satisfy all these inequalities, 
constituting a convex polytope. 

Any convex polytope admits two descriptions, one in 
terms of the half-space inequalities (H-description), and 
one in terms of its vertices (V-description) . The set Cd 
is, by definition, an H-polytope. For one qubit we have 
found the equivalent V-description: there are 6 vertices 
corresponding to the 6 basis vectors of the three MUB's 
used to define W. 

Let us now see how to find the equivalent V-description 
for the H-polytope Cd defined above. The starting point 
is the general expression for the d-dimensional Wigner 
function at phase-space point a (see il8j): 



Wo, 



(25) 



where the sum is over the probabilities associated with 
projectors corresponding to all lines Xij containing 
phase-space point a. By construction of the striations, 
each point belongs to exactly one line from each stria- 
tion. Thus for each point a, the sum above has {d + 1) 
terms pij, one from each MUB. 

We want to find which conditions on correspond 
to the demand of non-negativity of Wa for all points a 
and for all definitions of W using a fixed MUB set. For 



a fixed phase-space point a, changing the definition of 
W will correspond to picking a different set of (d -I- I) 
probabilities Pi J- in eq. (|25|l . one from each MUB i. There 
are only d^'^^^'^ ways of doing so, and this is the total 
number of expressions for Wa which we would like to 
take only non-negative values. These inequalities 
of the form Wa > constitute our H-description of the 
polytope Cd- 

From this H-description it is possible to find the equiv- 
alent V-dcscription using a convex hull program based 
on the QuickHuU algorithm While this is consid- 

ered to be a computationally hard problem in general, it 
was possible to do the calculation for d < 5. The results 
are similar to the one-qubit case. For d = 3,4 and 5 I 
found that the V-description of Cd has d(d -|- I) vertices, 
each corresponding to one of the MUB basis vectors used 
to define W. For d = 5, for example, the H-polytope 
C5 in 5^ — 1 = 24-dimensional p-space is delimited by 
5® — 15625 half-space inequalities corresponding to non- 
negativity of W in all definitions. The V-description of 
this polytope consists of exactly 5 x 6 = 30 vertices, 
each corresponding to one of the basis vectors of the six 
MUB's. 

It is not hard to see that the MUB basis vectors have 
non-negative W in every single definition of W using 
those basis vectors. For d < 5, the computation de- 
scribed above shows that these are the only pure states 
for which this is true for all definitions, and moreover 
that the mixed states with the same property are exactly 
their convex combinations. In the next section we will 
discuss the implications for quantum computing. 



A. The set Cd and quantum computing 



Given the results above for d 
ural to make a conjecture: 



2, 3, 4 and 5, it is nat- 



Conjecture: For any power-of-prime Hilbert space di- 
mension d, the polytope Cd is equivalent to the 
V-polytope whose vertices are the basis vectors of 
the MUB's used to define W. 

It is easy to show that the latter V-polytope is contained 
in Cd; the converse seems to be harder to prove for gen- 
eral power-of-prime d. 

To understand the relevance of this conjecture for 
quantum computation, we need to review some known 
facts about mutually unbiased bases. Various authors 
have come up with constructions of MUB's 30., 3&, 2^ 
ESEllil- For d = 2^, i.e. for the case of qubits, the 
corresponding basis vectors can always be written as si- 
multaneous eigenstates of tensor products of single-qubit 
Pauli operators More generally, for any power-of- 

prime dimension d the MUB basis vectors can be chosen 
to be eigenstates of generalized Pauli operators 39, 4^. 
In other words, the basis vectors of MUB's can always be 
chosen to be a set of stabilizer states. 
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The conjecture would then mean that for any d, the 
set Cd would be the convex hull of the particular subset 
of stabilizer states used as basis vectors of the MUB's. 
Pure-state quantum computation which is restricted to 
states in Cd would then only ever reach the MUB basis 
vectors, all of which admit an efficient classical descrip- 
tion using the stabilizer formalism. 

As I have already remarked, it is well-known that 
quantum computation that is restricted to pure stabilizer 
states and gates from the Clifford group is not universal, 
and can be efficiently simulated on a classical computer 
prf . If the conjecture holds, then a wide class of uni- 
tary dynamics restricted to pure states inside Cd would 
be easy to simulate on a classical computer - for exam- 
ple, any gate from the Clifford group. This suggests that 
non-trivial pure-state quantum computation may require 
states outside of Cd, i.e. having negative Wigner func- 
tion W in at least some of the definitions of W in the 
class proposed by Gibbons, Hoffman and Wootters [l^ . 

VI. CONCLUSION 

In this paper I related the negativity of the dis- 
crete Wigner function W with quantum computational 
speedup. I characterized the set of states with non- 



negative W using the class of Wigner functions W pro- 
posed by Wootters and collaborators in 0, 0| . Wigner 
functions in this class are defined using projectors on a 
fixed set of mutually unbiased bases (MUB's). The set of 
states with non-negative W depends on which definition 
of W from this class we choose. 

I defined the set Cd of d-dimensional states having non- 
negative W simultaneously in all definitions of W in the 
class. States in C2 are classical in the sense of failing 
to provide quantum computational advantage in a recent 
model proposed by Bravyi and Kitaev p^ . For dimen- 
sions 2 < d < 5 I showed that Cd is the convex hull of the 
stabilizer states used as basis vectors in the MUB set. 

These results for small dimensions d support the con- 
jecture that for any power-of-prime dimension d, the set 
Cd is the convex hull of a set of stabilizer states. This 
would mean that pure states in Cd are 'classical' in the 
sense of admitting an efficient classical description using 
the stabilizer formalism. This suggests that states with 
negative W may be necessary for exponential pure-state 
quantum computational speedup. 
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